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A major problem in exploiting microscopic systems for developing a new technology based on the 
principles of Quantum Information is the influence of noise which tends to work against the quantum 
features of such systems. It becomes then crucial to understand how noise affects the evolution 
of quantum circuits: several techniques have been proposed among which stochastic differential 
equations (SDEs) can represent a very convenient tool. We show how SDEs naturally map any 
Markovian noise into a linear operator, which we will call a noise gate, acting on the wave function 
describing the state of the circuit, and we will discuss some examples. We shall see that these gates 
can be manipulated like any standard quantum gate, thus simplifying in certain circumstances the 
task of computing the overall effect of the noise at each stage of the protocol. This approach yields 
equivalent results to those derived from the Lindblad equation; yet, as we show, it represents a handy 
and fast tool for performing computations, and moreover, it allows for fast numerical simulations 
and generalizations to non Markovian noise. In detail we review the depolarizing channel and the 
generalized amplitude damping channel in terms of this noise gate formalism and show how these 
techniques can be applied to any quantum circuit. 



PACS numbers: 03.65.Yz, 03.67.Hk, 03.67.Lx 



I. INTRODUCTION 



The development of both theoretical and experimen- 
tal research on quantum information is opening the way 
to a new technology based on the quantum properties 
of microscopic systems; its potentialities are extraordi- 
nary, but the possibility of actually exploiting quantum 
properties for building new physical devices is not yet 
completely clear. The main reason for this limitation, as 
is well known, is that quantum systems are highly sen- 
sitive to the influence of the surrounding environment 
which tends to destroy quantum coherence: it becomes 
important to analyze to what extent external influences 
disturb the time evolution of quantum systems. 

The most common techniques employed in studying 
the interaction of a quantum system with an environ- 
ment mainly rest on the master equation approach [1], 
the operator-sum representation method [2] , and stochas- 
tic unravelings of master equations in terms of random 
quantum jumps [3] or stochastic differential equations 
(SDEs) [4]. The SDE approach has gained increasing 
popularity in recent years, but attention has focused 
mainly on non linear SDEs which in general arc difficult 
to work with. In this paper we show how SDEs can be a 
flexible and handy mathematical tool for analyzing many 
physical situations analytically and numerically. The key 
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property of SDEs we will use is that among the differ- 
ent stochastic unravelings of a given master equation of 
the Lindblad type, there is always one that is linear [5]: 
by resorting to this specific unraveling, the power of the 
superposition principle can be used to analyze the evo- 
lution of the open system. As we shall see, the effect 
of the environment can then be described in terms of a 
linear and stochastic matrix, which can be manipulated 
like any other standard quantum gate, when any quan- 
tum protocol is analyzed. For this reason we will call it 
a noise gate. 

This approach has some advantages with respect to the 
other ones, at least in certain circumstances: 

1. During computations it allows one to work with 
state vectors instead of density matrices, even if 
the system is open; this makes the analysis simpler, 
since it allows the system to be treated as if it were 
closed, the effect of the environment being modeled 
by a random potential. 

2. It is predictively equivalent to the Lindblad ap- 
proach in the limit of Markovian interactions; at 
the same time, it can be generalized to non Marko- 
vian dynamics [6], which are a subject of increasing 
theoretical interest [1,7] for the description of sev- 
eral important physical phenomena [8]. 

3. In many important cases it allows exact mathemat- 
ical results to be computed; also those, such as the 
long-time behavior or the limit for a large number 
of qubits, which cannot be computed numerically. 
When exact results cannot be obtained, it allows for 
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a perturbation analysis or for fast numerical simu- 
lations [9]. Alternatively, in some cases, part of the 
analysis can be done analytically and part numeri- 
cally, thus simplifying the overall work. 

The paper is organized as follows. Section II sets up the 
general formalism which allows the effect of the environ- 
ment to be expressed through noise gates. In Sec. Ill we 
compute the corresponding noise gates for four standard 
noise channels. Sections IV- VI contain pedagogical ex- 
amples of how the noise gate formalism can be applied to 
any quantum circuit. In Sec. VII we conclude with final 
remarks. 



II. MASTER EQUATIONS, SDES AND NOISE 
GATES 



Let us consider a master equation of the Lindblad type: 

j t p{t) = -l [H,p(t)} +1 Lp(t)tf -l{L%p{t)}, (1) 

where the Lindblad operator L (for the sake of simplic- 
ity, we consider only one such operator) summarizes the 
effect of the environment on the quantum system, and 
7 is a coupling constant. It is well known that Eq. (1) 
allows for different unravelings in terms of SDEs; what 
perhaps is less known is that, among such unravelings, 
there is always one which is linear [5], namely, 



d\tp t 



Hdt + i^jLdW t 



1 



itfLdt 



lV>t>, (2) 



where W t is a Brownian motion defined on a probabil- 
ity space (fi,.F,P). When L is a self-adjoint operator, 
Eq. (2) preserves the norm of \ipt)> however this is no 
longer true for general operators. In such cases one can 
always replace the above equation with a non linear one, 
which is norm preserving [10]. However, this is not nec- 
essary since, even for the unnormalized state \ipt) being 
a solution of Eq. (2), one can easily prove that 



P (t) EE E[|Vt><V*|], 



(3) 



which means that, when the stochastic average E (with 
respect to the measure P) is computed, the predictions 
of Eq. (2) are equivalent to those of Eq. (1). In this sense 
Eq. (2) is an unraveling of Eq. (1). 

Now, because of linearity, the solution of Eq. (2) can 
be generally expressed as follows: 



(4) 



and the linear operator N(t,to), which from now on we 
will refer to as the noise gate, can be treated like any 
other quantum gate, except for the fact that in general 
it is not unitary. Such a gate of course depends on the 
noise W t , which means that for each realization of the 
noise the system evolves following different "histories" . 



So all one has to do is to solve Eq. (2) to get the correct 
expression for the noise gate N(t,to), which can then 
be inserted in a quantum circuit in the usual way; after 
having computed the square modulus of the probability 
amplitudes, their stochastic average can be calculated to 
obtain the correct physical predictions. 

In the following we give the exact expressions of four 
important quantum noise gates, and on the basis of two 
example quantum circuits we show how the proposed 
noise gate formalism can be used to treat the effects of 
quantum noise in any quantum circuit. 

A note before concluding. In the Introduction we men- 
tioned that other (infinitely many different) unravelings 
of Eq. (1) in terms of SDEs are possible [4]; in such 
cases Eq. (2) is replaced by other, structurally different, 
SDEs, which in general depend in a non linear way on 
\ipt)- Despite this, the relation (3) still holds true for any 
such equation, implying that, at the statistical level, the 
predictions computed through the SDEs are equivalent 
to those computed through the master equation. This 
also implies that, at the statistical level only, all effects 
of non linearity are " washed away" . The disadvantage of 
such unravelings with respect to ours is that, precisely 
due to their non linearity, the solution of the SDE can 
not be written as in (4) in terms of a linear operator, and 
therefore the superposition principle can not be used to 
infer the effect of the noise on a generic state, once its 
effect on a basis is known. 



III. EXAMPLES OF NOISE GATES 
A. Bit Flip, Phase Flip, Bit-Phase Flip Channels 

These three channels are accurately described in [2] 
within the framework of the quantum operator-sum rep- 
resentation. It is not difficult to represent them in terms 
of a master equation of the form (1): the corresponding 
Lindblad operator L turns out to be one of the Pauli 
matrices, a x for the bit flip channel, o~ z for the phase flip 
channel and a y for the bit-phase flip channel. 

For each of these operators, the SDE (2) becomes 



■AWt 



1 



-fldt 



\4n. 



(5) 



where I is the identity matrix and k = x,y,z. Since 
the matrices appearing in Eq. (5) obviously commute, it 
can be solved by means of standard techniques [12]; the 
solutions are: 

AWi(Mo) = exp[iV7MW t - W t0 )], (6) 
AWMo) = exp[i^ya z {W t - W to )] , (7) 
N mt . Phm (t,t ) = exp[i^ya y (W t -W to )]. (8) 

As we see, the above gates lead to a nice physical inter- 
pretation of the effect of the environment on a qubit: it 
randomly rotates the qubit along a specific direction, the 
randomization being proportional to the strength of the 
coupling constant 7. 
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B. Amplitude Damping Channel 



The noisy CNOT gate acts on the qubits as follows: 



The amplitude damping channel is also described in [2] , 
and the associated Lindblad operator is cr~ = |0)(1|; 
written in terms of the components a t = (0\ipt) an d 
Pt = {l\i/) t ), Eq. (2) becomes: 

da t = yfiPtdWt, (9) 
dfit = -( 7 /2)/3 t dt (10) 

The solution, expressed in matrix notation, is: 

^wmo) = a if (A y (ii) 



|0,0) 



(1) 



^Vl(*2,*l) ® JV2(*2,*l) |0,0) 



e 2 
■t 



^7 / e~i s dW a 

Jt 



(12) 



This channel models loss of energy to the environment: 
the |1) state decays to |0) at a given rate 7 while |0) is 
stable. 

We now turn to first applications of this formalism, 
during which we will spot some general features of the 
noise gates, which are useful for simplifying the calcula- 
tions. 



IV. APPLICATION 1: NOISY C-NOT GATE 

As a first example, we analyze the controllcd-NOT gate 
(CNOT) by assuming that, before and after its applica- 
tion, the involved pair of qubits are subject to environ- 
mental noise, as shown in Fig. 1. Such a quantum cir- 
cuit is interesting because the final state of the target 
qubit will depend not only on the noise acting on it but, 
through the CNOT gate, also on the noise acting on the 
control qubit: the noise gate formalism allows this depen- 
dence to analyzed in simple terms. For dcfinitcncss, let 
us take |0, 0) as the input state. Let us moreover assume 
that the noises acting on the two qubits are independent 
of each other; this assumption is of course justified only 
if, e.g., the two physical states encoding the two qubits 
are separated by more than the correlation length of the 
noise, so that the surrounding environment acts indepen- 
dently on them. 



*i>- 



*2)~ 



N-2 



Ni T N i - l*i 



N- 2 



- *2) 



time 
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FIG. 1: CNOT gate for two qubits subject to environmental 
noise. We assume that the two noises represented by the gates 
Ni and N2 are independent. 



(2) 



(3) 



J2nf' 0) (t 2 ,h)ni j ' \t 2 , tl )\i,j) 

i,j=0 
1 

£ n^ 0) (t 2 ,M4 J ' 0) (Wi)K,*©J> 

ij=0m,t0 

■n ( r\t 3 Mn ( 2 - m \t 3 ,t 2 )\ m ,£), (13) 

where n& (%,t a ) is the (i, j)-th coefficient of the matrix 
N a (tt, t a ), a — 1, 2. Step (1) takes into account the effect 
of the environment from time t\ to time t 2 , after which 
the CNOT gate is applied: the latter corresponds to step 
(2). In step (3) we assume that the environment acts on 
the qubits until time t 3 . Let us call |0, 0) the final state. 

For a closed system we would have of course |0, 0) = 
|0, 0), but with the environment interacting with the two 
qubits, 1 , 0) becomes a random entangled state; to un- 
derstand the effect of the noise, let us compute the fidelity 
F = E[|(0,0|070)| 2 ] of the noisy protocol (the stochastic 
average E takes into account all possible realizations of 
the noise). Because of the Markov property of Wiener 
processes, we have the following property: 

Property 1. Two gates N(t 2 ,ti) and N(t 4l t 3 ) are in- 
dependent whenever (ti,t 2 ) fl (t 3 ,t 4 ) is empty. This fol- 
lows directly from the fact that a Brownian motion has 
independent increments. 

According to this rule, and keeping in mind the assump- 
tion that the noise gates acting on the two qubits are also 
independent, F is a combination only of terms having 

the form F,[ria' k \tb, t a )na ' k \tb,t a )*] with a = 1,2 and 
(tb,t a ) = (t 2 ,ti) or (t 3 ,t 2 ), since all other terms vanish 
when averaged, because of the statistical independence. 
Now it is just a matter of choosing the type of noise that 
best describes the environment, and computing the re- 
quired expectation values. 

Let us consider, as an example, the bit flip gate given 
in Eq. (6): 

with AWf = Wt — Wt ; one easily verifies that, for a 
standard Wiener process, the following equalities hold 
true: 

E[cos 2 (^(W tb - W ta ))] = Pa (t b - t a ) (15) 

1 +cxp[-27 a (t b - t a )] 
2 

E[sin 2 (^(W tb - W t J)] = p a (t b - t a ) (16) 

= 1 -p a (t b - t a ), 
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while E[cos(V7a(W tl> - W t J) sm(^(W tb - W ta ))] = 
0. From these expressions one can derive the following 
expressions for the correlation functions of the coefficients 
of the noise gate: 

E[n£ k \t b) t a )n<Z'> k '){t b ,t a )*} = Pa (t b -t a ) (17) 

if i = k, i = k ; 

E[ng'*)(t 6 ,*a)nS" ,fc ' ) (t6,tB)*] =P«(*6-ta) (18) 

x£i^k,i'jt k'; 

E[n£ k \t b ,t a )ng' k '\t b ,t a y] = (19) 

in all other cases. 

Given this, one easily gets the following expression for 
the fidelity F as a function of the coupling constants 7 Q 
and of the time intervals during which the noises act on 
the CNOT gate: 



F = pi(t 2 


- h)p 2 {t 2 


- h)pi(t 3 


- h)p 2 {h 


-h) 


+ Piih 


- h)p 2 {h 


- h)Pi{h 


- h)p 2 {h 


-h) 


+ Pi(t 2 


~ h)P 2 {t 2 


- h)p\{h 


- t 2 )p 2 {h 


-t 2 ) 


+ Pi{h 


- ti)p 2 (h 


- h)Pi{h 


- t 2 )p 2 (t 3 


-*a) 



(20) 



In particular, if we assume that the two noises have the 
same strength (71 = 72) and that the time intervals dur- 
ing which they act are the same £3 — t 2 = t 2 — 1\ = T, we 
obtain the simplified formula: 

F(T) = Ap{tf - 5p(T) 2 + 2p(T), p(T) = i-tl , 

(21) 

which shows, as expected, that F starts from 1 and de- 
creases exponentially in time to 1/4: the formula displays 
the whole time evolution. 



V. APPLICATION 2: TRANSFER OF AN 
ENTANGLED STATE THROUGH A SPIN 
CHAIN. 

One of the most common problems in quantum infor- 
mation theory is the transfer of information through a 
noisy channel, which is often analyzed by modeling the 
channel with a spin chain [11]. By means of the spin 
chain (see Fig. 2) we want to demonstrate the power of 
the noise gate formalism for analyzing whole quantum 
circuits. For the spin chain we chose a standard model 
consisting of a chain of n + 1 qubits, such that the first 
two qubits are in a given normalized entangled state |?/>), 
while the remaining qubits are in a general normalized 
state \<p); then the global initial state is |^) = \rp) (g> \<f») 
with 

1 

IV^> = E ^o.uKon) (22) 

i o ,ii=0 
1 

\4>) = X/ C i2,...U*2 ■ • -in)- (23) 
i 2 ,...in=0 



The state of qubit 1 can be transferred to the far end 
of the chain by means of a sequence of n — 2 swap op- 
erations; for a circuit not subject to noise, at the end of 
the protocol the first and last qubits are in the entan- 
gled state \tp), factorized from the rest of the chain: we 
now analyze how the protocol is changed by the effect of 
a noisy environment, which we describe by n + 1 noise 
gates Nk(ti,tj), k — 1, . . . ,n, each acting on a different 
qubit for the time interval (tj,ti). We assume that the 
environment is random enough to act independently on 
each qubit; this means that we assume that the noise 
gates are independent. 

The time evolution of the global state can now be im- 
mediately computed: due to the linearity of the noise 
gates, and with reference to Fig. 2, one easily gets for 
the final state at time t = t n 

l*)= Yl ^oM^,...in\io)®\h)®---K)> (24) 

io,~ in 

where the states \ik) are defined as follows: 

\io) = N (t n ,t )\i ), 

|«fe) = N k (t n ,tk)N k+ i(tk,t ) \ik+i), fc = l,...n-l 

n-l 

K) = JJ iV„_ fc (t„_ fc ,t„_ fe _i) (25) 

fc=0 

so the problem is mathematically solved. In order to test 
the effect of the environment on the transmission pro- 
tocol, we compute the reduced density matrix referring 
to the 0th and nth qubits, obtained from the full den- 
sity matrix E[|\&)(\&|] by tracing away all other degrees 
of freedom (from 1 to n — 1): 

P(o,«)(t») = Tri,...,(n-i)E[|*)(*|] (26) 
= E[\i>)(i>\-Tr\&m, (27) 

where we have defined: 
1 

\i>) = ^0,^0(^0) |*o)iV|ii) 

io,ii=0 
1 

\4>) = E c i2 ,... in N x \i 2 )...N n _ l \i n ), (28) 
i 2 ,...in=Q 

and we have used the short-hand notation: 

N k = N k (t n ,t k )N k+1 (t k ,t ) (29) 

n-l 

N = JJjV n -fc(*n-fc,*n-fc-l). (30) 
fc=0 

Note that the state \-tjj) is a linear combination of the two 
states |«o) and \i n ) (see Fig. 2) — this is the reason why 
the trace does not affect — while \<fi) is a linear 

combination of the remaining states . . . \i n -i)- 

In our setting, Property 1 together with our assump- 
tion of independence of the environments lead to the in- 
dependence of the statistics of the noise gates Ni(t v ,t u ) 
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FIG. 2: Quantum circuit for the transmission of an entangled state through a spin chain in a noisy environment. 



and Nj(t s ,t r ) whenever i ^ j or (t u ,t v ) n (t r ,t s ) is 
empty. One easily verifies that No(t n ,to), N and iV fc 
for A: = 1, ... ,n — 1 are independent; this can be im- 
mediately checked in Fig. 2. This also means that the 
statistics of IVvW is independent of that of \(j>)(<p\ : and 
their average values can be computed separately: 

p (0 ,„)(i„) = E[|^}(^|] ■ E[Tr|0}(0|]. (31) 

Another important property of the noise gates, which 
we shall now use to simplify the above formula, is the 
following. 



Property 2. Given an initial normalized n-qubit state 
|V>o) which evolves, according to a SDE of the type (2), 
to a random state \ipt) — iV(t, io) l^o) j then the following 
equality holds true: TrE[|?/> t )(?/> t |] = 1 for any t. This 
property is a direct consequence of Eq. (3) and of the 
fact that pit) satisfies Eq (1), which is of the Lindblad 
type and thus trace preserving. 



For the sake of brevity, we denote fa, . . . ,i n by i; the full 
expression of TrE[| \<j>) {<j>\\ is, according to Eq. (28), 



TrE[| |0) (4>\] =X>[£°»4 II 0*l^fc-i(*«.**-i)^fc(**-i»*o)|i*> (AW* k (t k -ut )NLi(tn,tk-i)\jk)]; (32) 

j k=2 

we now insert two identities |Z) (l\ between the noise matrices, and after a rearrangement of the terms we get 



TrE[| \$) (4>\] = n TrE K-i(*»' tfc -i) ( l 'k\ N k-i{tn,t k -i)] 'EL-4 (l k \N k (t k _ u t )\i k ) (if k \N*(t k ^,t )\l' k 



iA'Ll' k=2^- 



(A) 



(33) 



where the factorization of the two average values is again 
justified by the assumption of independence of the envi- 
ronments and by Property 1. 
By Property 2, we have 



(A) = TvE[\l k ) {l' k \ 

so Eq. (33) simplifies as follows: 
TrE[| \$) (4>\] = 



(34) 



(35) 



TrE [E IT CiC%N k {t k -x,t Q ) \i k ) Ntftk-uh) 

k=2 



which, again by Property 2, gives 



TrE[| 



= Tr 



= 1. 



(36) 



Accordingly we are left with the expected simple result: 



P(0, n) (*n) = E[|$$|], 



(37) 



from which any relevant piece of information can be ob- 
tained. 
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A. Fidelity of the transmission protocol 

As an application of this formula we now compute the 
fidelity 



F = Tr[p (0:n) (i n )^)(V|]-E[|(V'l^)| 2 



(38) 



of the transmission protocol; here, again, the noise gates 
come in handy in the computation since we may work 
with random vectors instead of density matrices. In order 
to focus our attention only on the effect of the noise on 
the qubit that has been transmitted, we neglect the effect 
of the noise gate N (t n ,t ) on the Oth qubit. We denote 
the random matrix components of N by 



N 



a b 
c d 



(39) 



and compute 



= C 8 *)* O'oii 1 1 <8> N \i h) (40) 



Aa + Bb + B*c + (1 - A)d, 



(41) 



where 
A := 



x o,oi 



l i,ol ) 



B 



«o,o«o,i + <o°i,i- ( 42 ) 



In order to become more concrete we chose the ampli- 
tude damping gate for iVfc, cf. Eq. (12), 



Nk(tk,tk-i) 



, Cfc dk 



1 WTkJt 1 '-*'<™> 



v0 



e 2 



l) 



(43) 
(44) 



The coupling constant 7^ represents the strength of the 
interaction of the fe-th noise with the fc-th qubit. 

By definition of N, cf. Eq. (30), its matrix components 
are given by 



0. 



(45) 



fe=i 



while b = b n (for n qubits) is defined by the recursive 
formula b n := b n + d n b n -i, with b\ = b\ where b\ = 
b\. Now we have all we need to compute the fidelity of 
this protocol. Plugging these matrix components into 
Eq. (41) we get by Eq. (37) 

F = E[\Aa + Bb+(l- A)d\ 2 ] (46) 
= A 2 + \B\ 2 E[\b\ 2 ] +(1- Afd 2 . (47) 

In the last step we have used the fact that only b is ran- 
dom and that only the terms quadratic in b give a non 
zero contribution to the expectation value. Using the re- 
cursive definition of b, cf. Eq. (45), and again collecting 



only the terms quadratic in the random variables we 
compute 

E[|6| 2 ] =E[\b n + dJ n - 1 \ 2 ] (48) 

= E[\b n \ 2 ] + d 2 Mb n ^\ 2 ] (49) 

n 

= l-e- r , T = J27*(tcx-t a -i) (50) 

a=l 

— r 

by induction. Together with d = e~2 we arrive at the 
formula: 



-FAmDa ~ 



A+(1-A)e-* +|B| 2 (l-e- r ) (51) 



For the bit flip, the phase flip and the bit-phase flip 
channels, cf. (6)-(8), we use the noise gates 

Niitut^x) = cxp (i^yia K (Wt, - WVJ) (52) 

with k = x,z, y, respectively, and so compute N for the 
three cases according to equation (30): 



AWi = 



Ny, 



Ah 



,i<j> 



cos (p 
— sin < 



sin c 

COS( 



for 



E 

fc=0 



Using equation (40) the fidelity turns out to be 
F BWP , = E (cos 2 6 + 4ReB 2 sin 2 6) 



F PhM =E(|l + A(, 



,2i<t> 



1)| J 



-^BitF 



= E(cos 2 <i + 4Im J B 2 sin 2 . 



(53) 
(54) 
(55) 

(56) 



(57) 
(58) 
(59) 



Here Re and Im denote the real and imaginary parts 
respectively. After evaluation of the expectation value 
the fidelity is given by 



F K = 



1 - n 2 

1 9 K -2r 



(60) 



2 2 

where k denotes Bitpi, PhFi BitPhFi, and g mtF1 = 2ReB, 
gphPi = 2^4 — 1 and g B it-PhPi = 2ImS. As expected, 
the fidelity decreases exponentially in time, reaching an 
asymptotic finite value which depends both on the initial 
entangled state and on the type of noise. More gener- 
ally, the above formula displays the full dependence of F 
on the different parameters entering the protocol, in par- 
ticular on the time between two subsequent application 
of a swap operation and on the strength of the different 
noises. It then applies, e.g. to non homogeneous envi- 
ronments, where some qubits feel a stronger decoherence 
effect than others. One can easily generalize the above re- 
sult by including also uncertainties in the times at which 
the different swap operations are applied. 
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VI. LINEAR COMBINATIONS OF NOISES 

So far we have looked only at SDEs that are explic- 
itly solvable. In this section we want to consider more 
complicated noise channels for which an explicit solution 
might not be available. In the chosen examples of this 
section we shall see that this is already the case when we 
combine two or more of the noise channels that we have 
discussed so far. 

In general, a linear combination of noises acting on the 
same quantum system can be described by 



yW L K dWi K) - \^LlL K dt 



IV>*>, 



(61) 

where in contrast to (2) we have neglected the Hamil- 
tonian but in addition have multiple Brownian motions 
Wt . For Lindblad operators L K this in turn leads to 
the corresponding master equation 



dt 



p(t) = E 



l iK) L,P(t)Ll 



{L{L K ,p(t)} 



(62) 



in the sense that (3) holds accordingly. 

As an important feature of the noise gate formalism 
we notice that, since quantum averages are always ex- 
pressed as the square modulus of the scalar product of 
two vectors, the coefficients of the noise gates always en- 
ter the stochastic averages in quadratic combinations, 



which with a little abuse of terminology we shall refer 
to as second moments. Now, although (61) might not 
in general be explicitly solvable it is often possible to 
infer from it all second moments of the noise gate that 
it describes either analytically or numerically; see [12] 
for a discussion of this topic. As we shall demonstrate, 
this can be done in an easy way whenever a solution 
to the corresponding master equation is available. Hav- 
ing these second moments computed either analytically 
or numerically, one may still work in the noise gate pic- 
ture even without having the explicit form of the noise 
gate, which in many circumstances can be more intuitive 
and faster. For the following discussion let us denote the 
(i,j)-th unknown coefficient of a noise gate N(tb,t a ) by 
n 1 {tbit a ), i,j — 0,1. Then, if the second moments of 

this noise gate, i.e. E (tb,t a )n( kl ' ) (tb,t a )*^J for any 

k,l — 0,1, are available one may perform any com- 
putation of a quantum average using the noise gate for- 
malism and in the end plug in the second moments when 
evaluating the stochastic average. 

In order to compute the second moments whenever 
a solution of the master equation is available, consider 
\ip(t a )) to be the initial value of the SDE and p(t a ) — 
\ip(t a )) (V>(i a )| the initial value of the master equation, 
both at time t a . As discussed before, the solution to the 
SDE at time tb can be expressed via the noise gate it de- 
scribes as N(tb,t a ) \i/)(t a )). By (3) the two entities p(tb) 
and E(N(t b ,t a ) \ip(t a )) {ip(t a )\N(t b ,t a )*) must be equal. 



For \lp(t a )) = J2i a i K) we nave 

Poo(tfe) = M 2 E (|4° 0) | 2 ) + l«i| 2E (Vt 1 '! 2 ) + 2Rea aiE (nf^n^*) (63) 

p Q1 (t b ) = |a | 2 E (n£ 0) n£ 0) *) + K| 2 E (n^n^*) + a a*E (n^n^*) + a^E (n^rc^*) = Pl0 (t b )* (64) 

p n (t b ) = \a \ 2 E (|r4! 0) | 2 ) + | ai | 2 E (|4 U) | 2 ) + 2Rea aiE (n$ 0) n^*) , (65) 

where Pij(t) = (i\p(t)\j) and T = tb — t a . Coefficient comparison then easily leads to the second moments. 



In the following we apply this scheme to the spin chain 
of the previous section treating two prominent represen- 
tatives of combined noise channels which are known as 
depolarizing and generalized amplitude damping channel, 
see [2]. 



terms of Eq. (61), k — 1,2,3 and the L K are the Pauli 
matrices a x , <r y , <r z . Here the effect of the environment is 
to randomly rotate the qubit around the x, y, z axis, with 
the randomization being proportional to the strength of 
the coupling constants 7^, 7^. 



A. Depolarizing Channel 



The depolarizing channel is a linear combination of 
the bit flip, phase flip and bit-phase flip channels. In 
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Hence its master equation takes the following form: 



d . . 
for which 



/ ( Pll W-pooW)7 (1 ' 2) 
UiW(7 (1) -7 (2) )-Pio(t)(7 (1) +7 ( 



PioW(7 (1) - 7 (2) ) - Poi(t)(7 (1) + 7 (2) + 2 7 (3) )\ 
7 (2) +2 7 (3) ) (p oW-Pii(i))7 (1 ' 2) )' 



Poo(M = i (|a | 2 (l + e- 2 ^ (1 * 2) ) + | ai | 2 (l - e- 2 ^ (1 ' 2 ')) 

= ^ (aoaUe- 2 ^' 1 ' 3 ' + e^ T ^) + ala^^ - e^^)) = p 10 (t b )* 
Pn(h) = \ (l«o| 2 (l - e- 2 ^ (1 ' 2> ) + | ai | 2 (l + e- 2 ^ 1 ' 2 ')) 

is the solution for initial value p(t a ), where we have used T = t b — t a and ^ m ' n "> = + ~f( n \ 

I 



(66) 

(67) 
(68) 
(69) 



By coefficient comparison with Eqs. (63), (64) and (65) 
one finds 



r la 

Hi 



-2T7 ( 
-2T7 ( 



1 j- e -2T7 (2 ' 3) 



k = i' 
k^i' 



-2T7 



(2,3) 



e -2T7< 1 ' 3 ») ti) k,i',k'= 0,0,1,1 
e -2T 7 (1 ' 3) ) t i^,i',k' = 0,1, 1,0 



. else 



(70) 



defined in the previous sections and the inverse process 
for which |1) is stable and |0) decays. In terms of Eq. 
(61), K = 1, 2 and the L K are the operators cr_ and a + = 
|1)(0|. The amplitude damping channel that we have 
discussed in the previous sections is the zero temperature 
limit of this channel. For non zero temperature the qubit 
may now also gain energy at the rate 



In order to apply this noise channel to the spin chain 
circuit discussed above we only need to plug these terms 
into Eq. (38) where we use the same abbreviations as 
in Eq. (40). We shall label the coupling coefficients for 

the a-th qubit by *yffl , , 7a? ^ for all a = 1, . . . , n when 
evaluating the product in (30). The computation of the 
fidelity is then straight forward: 



^P =2( A2 + ( 1 - A ) 2 )( 1 + e_2r(1 ' 2) )+ 
+ A(l-A)(e^ 3) +e^ 3) )+ 



+ \B\ 2 (l-< 



- 2r<1 ' 2,) )+Rei3 2 (e- 2r<2 ' 3)) ~ e - 2r(1 ' 3)) ) 



(71) 



where = £^ =1 ( 7 £ m) + 7^)^ - t a -i) and 7 £> 

denotes the ith coupling constant of the ath noise gate 
in the circuit. Note that the formula reduces to F K , Eq. 
(60), when all the coupling constants are set to zero ex- 
cept the ones associated with one Pauli matrix. Figure 
3 displays the time evolution of the fidelity of the spin 
chain under the influence of the depolarizing channel for 
a specific class of initial states. 



B. Generalized Amplitude Damping Channel 

The generalized amplitude damping channel in turn is a 
linear combination of the amplitude damping channel as 




time 



100 



FIG. 3: Time dependence of the fidelity F of the qubit pair 
\/A|01) + s/i — A 10} separated by the described spin chain 
circuit (see figure (2)) consisting of 100 qubits under the in- 
fluence of depolarizing channels, see Eq. (71). The time in- 
tervals between the swap operations were chosen to be equal 
to 1 while the coupling coefficients 7^ ' were chosen according 
to a Gaussian distribution centered around the 50-th qubit 
with covariance cov^ such that covi : COV2 : COV3 = 1:2:3. 
Due to the chosen class of initial states, the curvature of the 

surface along the time axis is due to the coupling of 7„ and 

(2) (3) 

7a while 7„ determines the curvature along the A axis. 
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This time its master equation takes the following form: 



dt 



_Po^) (7 (l) +7 ( 2)) ^ 

- £ #(7 (1) +7 (2) ) PooW7 (2) -Pii^7 (1) , 



d p(t)= (Pn( f )7 (1) -pooW7 (2) 



(72) 



for which 



ia r ^ (1)+ ^ rT + | ai r ^ (1) v rr) 

* -ET 

a5aie 2 - 1 



* -It 
a aje 2 - 1 



(73) 



is the solution for initial value p(t a ), where we used T = % — t a and T = 7W + 7^). 

I 



Again, by coefficient comparison with Eq. (63), (64) and 
(65) one finds 



E(M°°)(i 6 ,ia)| 2 ) = 


7 (i) 


+ 7 (2) e -rT 


(74) 




r 


E(|7i( 01 )(i b ,; a )| 2 ) = 


7 (i) 

T~ 


e- TT 


(75) 


E(\n^ 10 \t b ,t a )\ 2 ) = 


7 (2) 

T~ 


er TT 


(76) 


E(\n^\t b ,t a )\ 2 ) = 


7 (2) 


+ 7 (i) e -rr 


(77) 




r 



and 



E(n^°\t b ,t a )n^(t b ,t a ) 



(78) 



while all other second moments are equal to zero. As 
we have done with the depolarizing channel we apply 
this noise to the spin chain circuit discussed above and 
therefore we, again, only need to plug these terms into 
Eq. (38) using the same abbreviations as in Eq. (40). 
In order to keep the displayed formulas short we choose 
the coupling constants to be the same for all qubits, i.e. 
7a = 7W and 7i 2 ^ = 7 1 - 2 ' for all a — l,...,n, when 
computing the product in (30). We then find 



(A 



2 7 



(i) 



+ (1-^) 



2 7 



(2) 



(A 2l — + {l-Af< Y 
2A{1 - A)e^ (t "-* o) 



r 

B\ 2 )e- r ^ 



-to) 



(79) 



for r = ^ 



7^ ; + 7^ 2 ' and 7 W = 7^, such that all noise 
gates N a in the circuit have the same coupling constants. 
Note that also this formula reduces to (51) if j( 2 > is set to 
zero. An example of the time evolution of the fidelity for 
a specific class of initial conditions under the influence of 
generalized amplitude damping is shown in figure 4. 



VII. CONCLUSION 

We have suggested the noise gate formalism as a handy 
approach for analyzing the effect of the environment on 



quantum algorithms; it is very intuitive as it allows the 
influence of the environment to be treated in terms of 
noise gates, which can be manipulated like any other 
quantum gate. In many situations it makes the compu- 
tation easier, either analytically or numerically. We em- 
phasize again that it is especially interesting for numeri- 
cal simulations because linear SDEs can be integrated by 
standard methods [9] . In contrast to solving the Lindblad 
equation numerically, which roughly scales quadratically 
with the number of degrees of freedom, the numerical in- 
tegration of (2) scales only linearly, even if the noises are 
dependent. Finally note that SDEs can also be general- 
ized to model non Markovian quantum noise. 

Acknowledgements. The work of A.B. has been 
partly supported by DFG (Germany). The work of 
D.-A.D. has been supported by the EU Grant No. 
ERG 044941-STOCH-EQ. 




FIG. 4: Time dependence of the fidelity F of the qubit pair 
\/A|01) + y/l — A 10) separated by the described spin chain 
circuit (see figure (2)) consisting of 100 qubits under the in- 
fluence of generalized amplitude damping channels, see Eq. 
(79). The time intervals between the swap operations were 



chosen to be equal 1 while 7 



(i) 



7 



( 7 «+ 7 C a >) 



= 3:1:4. 



The latter is reflected in the curvature along the A axis. 
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